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For phase transitions in dissipative quantum impurity models, the existence of a quantum-to- 
classical correspondence has been discussed extensively. We introduce a variational matrix product 
state approach involving an optimized boson basis, rendering possible high-accuracy numerical stud- 
ies across the entire phase diagram. For the sub-ohmic spin-boson model with a power-law bath 
spectrum oc w", we confirm classical mean-field behavior for s < 1/2, correcting earlier numerical 
renormalization-group results. We also provide the first results for an XY-symmetric model of a spin 
coupled to two competing bosonic baths, where we find a rich phase diagram, including both critical 
and strong-coupling phases for s < 1, different from that of classical spin chains. This illustrates 
that symmetries are decisive for whether or not a quantum-to-classical correspondence exists. 



Quantum spins in a bosonic environment are model 
systems in diverse areas of physics, ranging from dissi- 
pative quantum mechanics to impurities in magnets and 
biological systems [1]. In this paper we consider the spin- 
boson model and a generalization thereof to two baths, 
described by Hsh = —h ■ (t/2 + Hhnth, with 
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The two-level system (or quantum spin, with (Jx,y.z being 
the vector of Pauli matrices) is coupled both to an exter- 
nal field h and, via ax and (Ty, to two independent bosonic 
baths, whose spectral densities Ji(a;) = tt \^^^5{uJ—^;Jq) 
are assumed to be of power-law form: 

Ji[(jj) — 2tt ai Lu]r ^ u)' 



< < = 1 
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Such models are governed by the competition between 
the local field, which tends to point the spin in the h 
direction, and the dissipative effects of the bosonic baths. 
Indeed, the standard one-bath spin-boson model 



(SBMl), obtained for ay 



0, exhibits an in- 



teresting and much-studied [1-7] quantum phase tran- 
sition (QPT) from a delocalized to a localized phase, 
with {ax) = or 7^ 0, respectively, as ax is increased 
past a critical coupling ax,c- According to statistical- 
mechanics arguments, this transition is in the same uni- 
versality class as the thermal phase transition of the 
one-dimensional (ID) Ising model with l/r^"*"* interac- 
tions. This quantum-to-classical correspondence (QCC) 
predicts mean-field exponents for s < 1/2, where the 
Ising model is above its upper-critical dimension [8, 9]. 

Checking this prediction numerically turned out to 
be challenging. Numerical renormalization-group (NRG) 
studies of SBMl yielded non-mean-field exponents for 
s < 1/2 [4], thereby seemingly negating the validity of 
the QCC. However, the authors of Ref. 4 subsequently 
concluded [10] that those results were not reliable, due 
to two inherent limitations of NRG, which they termed 
(i) Hilbert-space truncation and (ii) mass flow. Problem 



(i) causes errors for critical exponents that characterize 
the flow into the localized phase at zero temperature, 
since (a^) 7^ induces shifts in the bosonic displacements 
Xq — {B^ -\-W^)/y/2 of the bath oscillators which diverge 
in the low-energy limit for s < 1 and hence cannot be ad- 
equately described in the truncated boson Hilbert space 
used by NRG [11]. Problem (ii) arises for non-zero tem- 
peratures, due to NRG's neglect of low-lying bath modes 
with energy smaller than temperature [12]. In contrast 
to NRG, two recent numerical studies of SBMl, using 
Monte Carlo methods [6] or a sparse polynomial basis 
[5], found mean- field exponents in agreement with the 
QCC. Nevertheless, other recent works continue to advo- 
cate the failure of the QCC [13]. 

The purpose of this Letter is twofold. First, we show 
how the problem (i) of Hilbert-space truncation can be 
controlled systematically by using a variational matrix 
product state (VMPS) approach formulated on a Wilson 
chain. The key idea is to variationally construct an opti- 
mized boson basis (OBB) that captures the bosonic shifts 
induced by {ax) 7^ 0. The VMPS results confirm the 
predictions of the QCC for the QPT of SBMl at T = 0. 
(Problem (ii) is beyond the scope of this work.) Second, 
we use the VMPS approach to study an XY-symmetric 
version of the two-bath spin-boson model (SBM2), with 
Cix — oiy. This model arises, e.g., in the contexts of im- 
purities in quantum magnets [14, 15] and of noisy qubits 
[14, 16], and displays the phenomenon of "frustration of 
decoherence" [14]: the two baths compete (rather than 
cooperate), each tending to localize a different compo- 
nent of the spin. As a result, a non-trivial intermediate- 
coupling (i.e. critical) phase has been proposed to emerge 
for s < 1 [15], which has no classical analogue. To date, 
the existence of this phase could only be established in an 
expansion in (1 — s), and no numerical results are avail- 
able. Here we numerically investigate the phase diagram, 
and, surprisingly, find that the perturbative predictions 
are valid for a small range of s and a only. We conclu- 
sively demonstrate the absence of a QCC for this model. 

Wilson chain. — Following Refs. 3, 11, which adapted 
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Wilson's NRG to a bosonic bath, we discretize the latter 
using a logarithmic grid of frequencies LOki oc A~'' (with 
A > 1 and k a positive integer) and map Hhath onto a 
so-called Wilson chain of (L — 1) bosonic sites: 
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Here n^i = b\ibf.i, with eigenvalue n^i, counts the bosons 
of type i on chain site fc; the detailed form of the hopping 
parameters tki, on-site energies eki (both cx A^*^), and 
the coupling rji between spin component ai and site 1, 
are obtained following Refs. 17, 18. To render a numer- 
ical treatment feasible, the infinite-dimensional bosonic 
Hilbert space at each site k is truncated by restricting the 
boson number to < nj,i < d}. (dj, < 14 in Refs. 3, 11). 
The standard NRG strategy for finding the ground 

state of ^sb'' — ' ^1'^ + ^bath^"* iteratively diago- 
nalize it one site at a time, keeping only the lowest-lying 
D energy eigenstates at each iteration. This yields a L- 
site matrix-product state (MPS) [19-21] of the following 
form (depicted in Fig. 1, dashed boxes): 

^ ^A"MAiM...^^-iK_i]|a)|n). (4) 

Here |ct) = |t), \ \) are eigenstates of CTa;; the states \n) — 
|rti, . . . , riL-i) form a basis of boson-number eigenstates 
within the truncated Fock space, with nki\n) = nki\n) 
and < riki < d^. For SBM2, Uk — {nkx,nky) labels 
the states of supersite k representing both chains. Each 
^'^[rife] is a matrix (not necessarily square, but of maxi- 
mal dimension D x D, with A'^ a row matrix and A^~^ 
a column matrix), with matrix elements (A'^[rifc])^^. 

The need for Hilbert-space truncation with small d^ 
prevents NRG from accurately representing the shifts in 
the displacements x^i — (bj^^ ~\- 5],^)/-\/2 that occur in the 
localized phase. This problem can be avoided, in prin- 
ciple, by using an OBB, chosen such that it optimally 
represents the quantum fluctuations of shifted oscillators, 
^'ki ~ ^ki ~ {xki)- While attempts to accommodate this 
strategy within standard NRG were unsuccessful [11], it 
was shown to work well [5] using an alternative represen- 
tation of SBMl using a sparse polynomial basis. 

VMPS method. — We now show that an OBB can also 
be constructed on a Wilson chain. To this end, view the 
state \G) of Eq. (4) as a MPS ansatz for the ground state 
of "H^^', that is to be optimized variationally using stan- 
dard MPS methods [19-21]. To allow the possibility of 
large bosonic shifts, we represent the A-matrix elements 
as [22-24] (Fig. 1, sohd hues) 
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FIG. 1: Depiction of the MPS Eq. (4), with each yl-matrix 
expressed in an optimal boson basis via A = AV [Eq. (5)]. 



Here in effect implements a transformation to a new 
boson basis on site fc, the OBB, of the form \nk) — 

En'Jo Vln, Wk) with < nfe < dopt. (For SBM2, is 
a rank-3 tensor.) This ansatz has the advantage that the 
size of the OBB, dopt, can be chosen to be much smaller 
(dopt ^ 50) than dk- Following standard DMRG strategy, 
we optimize the A'' and matrices one site at a time 
through a series of variational sweeps through the Wilson 
chain. As further possible improvement before optimiz- 
ing a given site, the requisite boson shift can be imple- 
mented by hand in the Hamiltonian itself: we first deter- 
mine the "current" value of the bosonic shift (xki) using 
the current variational state |G), then use it as starting 
point to variationally optimize a new |G") with respect to 
the shifted Hamiltonian ^ (6^ , b^ ) = H^^^ (6^^ , b'^^ ) , 
with b'f. - = bki — {xki)/V2. The shifted OBB protocol, de- 
scribed in detail in Ref. 18, allows shifts that would have 



required d'^ « 10^° states in the original boson basis to 
be treated using rather small dk (we used dk = 100). 

Spin-boson model. — We applied the VMPS method to 
SBMl (ay — hy = 0), with dissipation strength a = 
and fixed transverse field hz = 0.1, at T = 0. We focussed 
on the QPT between the delocalized and localized phases 
in the subohmic case, s < 1. Here, the controversy [4- 
6, 10, 13] concerns the order-parameter exponents /3 and 
(5, defined via {(Tx) oc (a — ac)*^ at h^ = and (ax) cx /i^* 
at a = ac, respectively. QCC predicts mean-field values 
/3mf = 1/2, Smf = 3 for s < 1/2 [8], whereas initial NRG 
results [4] showed s-dependent non-mean-field exponents. 

In Fig. 2a, we show sample VMPS results for (a^) vs. 
(a — Qfc) for s = 0.3 at h^ — 0, where ac was tuned to 
yield the best straight line on a log-log plot. The results 
display power-law behavior over more than 3 decades, 
with an exponent /3 = 0.50 ± 0.03. Deviations at small 
(a— ac) can be attributed to a combination of finite chain 
length and numerical errors of VMPS. Fig. 2b shows {(Jx) 
vs. at a — ac, and a power-law fit over 6 decades 
results in (5 = 2.9 ±0.2. Power laws of similar quality can 
be obtained for aU s > 0.2 [18, 25] (see [18], Fig. S7). 

The exponents /3 and S obtained from such fits are 
summarized in Figs. 2c, d. For s < 1/2 they are con- 
sistent with the mean- field values predicted by QCC, 
also found in Monte-Carlo [6] and exact-diagonalization 
studies [5], but are at variance with the NRG data of 
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S=0.3, h =0.1 , a =0.0346142, A=2, L=60, D=40, d =12 
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FIG. 2: VMPS results for the order parameter of SBMl near 
criticality. a) {a^} vs. (a — Qc) at h^; — 0, and b) {a^} vs. at 
a = Qc, on linear plots (insets) or log-log plots (main panels). 
Dashed lines are power-law fits in the ranges between the 
vertical marks. c,d) Comparison of the exponents and 5 for 
different s obtained from VMPS, NRG [4], mean-field theory, 
and, in d), the exact hyperscaling result 5 = (1 -|- s)/(l — s) 
which apphes for s > 1/2. (See also [18], Fig. S7). 



Ref. 4. Since both NRG and VMPS handle the same mi- 
croscopic model hI^^ defined on the Wilson chain, but 
VMPS can deal with much larger d'j^ values (< 10^° 
in Fig. 2) than NRG, the incorrect NRG results must 
originate from Hilbert-space truncation, as anticipated 
in Ref. 10. Indeed, artificially restricting dk to small val- 
ues in VMPS reproduces the incorrect NRG exponents 
(see [18], Fig. S6). 

Two-bath model. — We now turn to SBM2, a general- 
ization of the spin-boson model. Here, the two baths may 
represent distinct noise sources [14, 16] or XY-symmetric 
magnetic fluctuations [14, 15, 26]. Perturbation theory 
shows that the two baths compete: A straightforward 
expansion around the free-spin fixed point (a = h = 0) 
results in the following one-loop renormalization-group 
(RG) equations aX h — 0: 

I3{ax) = {l-s)ax-axay, l3{ay) = {l-s)ay-aj;ay. (6) 

For a = CKj. = tty, these equations predict a stable inter- 
mediate-coupling fixed point at a* = 1 — s, describing 
a critical phase. It is characterized by (a) = 0, a non- 
linear response of (tr) to an applied field h, and a finite 
ground-state entropy smaller than In 2, all corresponding 
to a fluctuating fractional spin [15, 27]. This phase is 
unstable w.r.t. finite bath asymmetry {ax ay) and 
finite field. It had been assumed [15] that this critical 
phase exists for all < s < 1 and is reached for any a. 

We have extensively studied SBM2 using VMPS; the 
results are summarized in the h = phase diagram in 
Fig. 3a and the fiow diagrams in Fig. 4. Most impor- 
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FIG. 3: a) Phase diagram of SBM2 in the s-a plane for 
h = 0, with dissipation strenght a = ~ cty. The critical 
phase only exists for s* < s < 1, and its boundary Oc — >■ oo 
for s — >■ 1^. (Ref. 18 describes the determination of the phase 
boundary and gives a 3D sketch of the s-a-h^ phase diagram, 
see Fig. S8.) b) Tranverse-field response of SBM2, (cr^) oc 



h 



i/s' 



for four choices of s and a, showing free (diamonds), 
critical (squares) and localized (triangles, circles) behavior. 

tantly, we find that the critical phase (CR) indeed exists, 
but only for s* < s < 1, with a universal s* = 0.75±0.01. 
Even in this s range, the critical phase is left once a 
is increased beyond a critical value a^s), which marks 
the location of a continuous QPT into a localized phase 
(L) with spontaneously broken XY symmetry and finite 
{(Tx,y)- This locahzed phase exists down to s = 0, Fig. 3a. 
It can be destabilized by applying a transverse field hz 
beyond a critical value h%{a)., marking the location of a 
continuous QPT into a delocalized phase (D) with unique 
ground state (see Ref. [18], Fig. S9). Finally, for s > 1 
we only find weak-coupling behavior, i.e., the impurity 
behaves as a free spin (F). 

In Fig. 3b (and Ref. [18], Fig. SIO) we show results 
for the transverse-field response, ((T^) cx h^J^ , which can 
be used to characterize the different zero-field phases. 
{az) is linear in hz in L {5' = 1), sublinear in CR {6' > 
1), and extrapolates to a finite value in F. For CR, a 
perturbative calculation gives 1/6' = (1 — s) + ©([1 — s]^) 



(b) s<s* 




CR QCl 



FIG. 4: Schematic RG flow for SBM2 in the a-h^ plane 
{hx — hy — 0). The thick lines correspond to continuous 
QPT; the full (open) circles are stable (unstable) fixed points, 
for labels see text, a) s* < s < 1: CR is reached for small a 
and hz = 0, it is separated from L by a QPT controlled by 
the multicritical QCl fixed point. Eq. (6) implies that CR is 
located at a* = 1 - s -I- 0[(1 - s)^]. For finite h^, a QPT 
between D and L occurs, controlled by QC2. b) < s < 
s*: both CR and QCl have disappeared, such that the only 
transition is between D and L. 
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[15] (confirmed numerically in Ref. [18], Fig. Sllb), while 
the linear response in L corresponds to that of an ordered 
XY magnet to a field perpendicular to the easy plane. 

From the VMPS results, we can schematically con- 
struct the RG flow. Fig. 4. There are three stable RG 
fixed points for s* < s < I, corresponding to the L, D, 
and CR phases. From this we deduce the existence of two 
unstable critical fixed points, QCl and QC2, controlling 
the QPTs (Fig. 4a). Eq. (6) predicts that, as s 1", 
CR merges with F; this is consistent with our results for 
S' which indicate (5' — )• oo as s — )■ 1~ (Ref. 18, Fig. Sllb). 
The behavior of the phase boundary ac in Fig. 3a sug- 
gests that QCl moves towards a = oo for s -> 1^. Thus, 
for s > 1 only F is stable on the h — axis. Conversely, 
from Eq. (6) and Fig. 3a we extract that, upon lowering 
s, CR (QCl) moves to larger (smaller) a. From the ab- 
sence of CR for small s we then conclude that CR and 
QCl merge and disappear as s — > s*^. Consequently, for 
s < s* we have only D and L as stable phases, separated 
by a transition controlled by QC2, Fig. 4b. The merger of 
CR and QCl at s = s* also implies that the phase bound- 
ary between CR and L in Fig. 3b at s* is vertical at small 
a (Ref. [18], Sec. V.C), because the merging point on the 
a axis defines the finite value of ac at s s*^. 

Taken together, the physics of SBM2 is much richer 
than that of a classical XY-symmetric spin chain with 
long-range interactions, which only shows a single ther- 
mal phase transition [28]. Given this apparent failure of 
QCC for SBM2, it is useful to recall the arguments for 
QCC for SBMl: A Feynman path integral representation 
of Eq. (1), with non-zero /i^, can be written down using 
eigenstates of both ax and Uz . Integrating out the bath 
generates a long-range (in time) interaction for a^- Sub- 
sequently, the az degrees of freedom can be integrated 
out as well, leaving a model formulated in ax only. Re- 
interpreting the ax values for the individual time slices 
in terms of Ising spins, one arrives at a ID Ising chain 
with both short-range and l/r^^'^ interactions, with the 
thermodynamic limit corresponding to the T — > limit of 
the quantum model. Repeating this procedure for SBM2 
with h ~ 0, one obtains a Feynman path integral in terms 
of eigenstates of ax and ay . Importantly, both experience 
long-range interactions and hence neither can be inte- 
grated out. This leads to a representation in terms of two 
coupled Ising chains. However, upon re-exponentiating 
the matrix elements, the coupling between the two chains 
turns out to be imaginary, such that a classical interpre- 
tation is not possible [29]. In other words, a Feynman 
path-integral representation of SBM2 leads to negative 
Boltzmann weights, i.e., a sign problem. 

Conclusion. — Our implementation of OBB-VMPS 
on the Wilson chain brings the Hilbert-space truncation 
problem of the bosonic NRG under control and allows for 
efficient ground-state computations of bosonic impurity 
models. We have used this to verify the QCC in SBMl 
and to determine the phase diagram of SBM2, which is 



shown to violate QCC. This underlines that symmetries 
are decisive for whether or not a QCC exists. A detailed 
study of the QPTs of SBM2 is left for future work. 

The results for SBM2 also show that the predictions 
from weak-coupling RG are not valid for all parameters 
and bath exponents, in contrast to expectations. This 
implies that studying a three-hath version of the spin- 
boson model, which is related to the physics of a mag- 
netic impurity in a quantum-critical magnet [15, 27], is 
an interesting future subject. 
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I. DISCRETIZATION AND MAPPING TO THE 
WILSON CHAIN 

The spin-boson model represents a prototypical 
quantum-impurity setup, with the bath consisting of non- 
interacting particles. As such it is amenable to the con- 
cept of energy scale separation present in the NRG^"'^. 
For this, the quantum impurity Hamiltonian of the spin- 
boson model is mapped onto a so-called Wilson chain, 
which includes two steps: (i) coarse graining of the bath 
(logarithmic discretization), followed (ii) by a mapping 
onto a semi-infinite bosonic chain with the spin-impurity 
connected to its starting point. 

The bath spectral function Ji{uj) of each bosonic bath 
i is assumed to be non-zero in the interval lo G [0,ajc], 
with ujc an upper cutoff frequency. The bath spectral 
function Ji{uj) is discretized then in energy-space into 



intervals [w- 



] , marked by the decreasing sequence 



(m = 0, 1, . . .) with Wo 



and lim„ 



0. 



Assuming two identical baths i e {x, y} that couple to 
the Pauli matrices Ui of the impurity, respectively, the 
discretized Hamiltonian has the form 



i—x,y m— 



2^ 



(SI) 

Here (Bmi) is the creation (annihilation) operator, 
respectively, of a free boson with energy , that is cou- 
pled to the impurity spin with strength 7^. Moreover, 




1/2 



J{x)dx 



(S2) 



We prefer a logarithmic discretization scheme over a lin- 
ear or power-law discretization, since the study of critical 
behavior requires very small energy scales to be resolved. 
This would require too many chain sites for linear and 
even power-law discretization schemes. Moreover, loga- 
rithmic discretization is ideally suited to represent scale- 
invariant physics near a quantum phase transition, and it 
has the advantage that characteristic NRG information, 
such as the energy flow diagram used to analyze the fixed 
points of the system, can be extracted from our VMPS 
results, if desired. (This will be elaborated upon in a 
separate publication''.) 



In this paper we use the improved logarithmic dis- 
cretization recently proposed by Zitko and Pruschke^ . As 
this achieves a more consistent description of the bath, 
it reduces discretization effects and hence allows to de- 
termine phase boundaries such as the critical coupling 
strength Uc more accurately. Thus we choose the dis- 
cretization intervals as 



c^cAi- 



(to = 1,2,3,...) 



(S3) 



with A > 1 Wilson's logarithmic discretization 
parameter^, and z G ]0, 1] an arbitrary shift^. By solving 
the differential equation in App. C of Ref. 5 analyti- 
cally, we obtain the following explicit expressions for the 
parameters in Eq. (S2): 
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Sm 



l-A" 



(l-(-s)lnA 



, -(1 + S)(,T1 + ^) 



1+s 



(l + .s)lnA 
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1+s 



(m = 1,2,3,...) 



(S4) 



f^P|(l - A-^(i+«)), 



27ra 
1+s 



(Ai+^ - l)A-('"+^)(i+^) (m = 1,2,3,.. .) 

(S5) 



Having discretized the Hamiltonian, the mapping onto 
the Wilson chain is done numerically using standard 
Lanzcos tridiagonalization. For the calculations in this 
paper, we use z = 1. 

We find that, with the discretization scheme of Zitko 
and Pruschke,^ the A-dependence of ac is much weaker 
than for the traditional discretization scheme''"^, i.e., Uc 
converges rapidly as A is decreased towards 1. Similar to 
standard NRG, critical exponents do not dependent on 
A, as shown in Fig. SI for the exponent /?. 



II. OBB-VMPS OPTIMIZATION PROCEDURE 

As discussed in the main paper, we use a MPS of the 
following general form: 

|G>= E E^°M^'M---^'^-M^L-i]k)|n). (S6) 
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FIG. SI: (cTj;) as function of (q — Oc) for SBMl, using three 
different choices of the discretization parameter A (and, cor- 
respondingly, different chain lengths L). The data illustrate 
that the critical exponent /3, obtained from power-law fits to 
this data, is essentially independent of A. The dashed line 
represents a power law {a.j^ oc (a — with mean-field ex- 
ponent /3mf ~ 1/2. 




FIG. S2: Update procedure at one site when sweeping from 
left to right. The matrices surrounded by the dashed lines 
are the outcome of singular value decomposition (SVD). The 
filled squares or triangles indicate the "current focus" of each 
step. Filled arrows indicate the local update loop (highlighted 
by the shaded area), open arrows indicate its entry and exit. 



The A-matrix elements are represented as 

fifc=0 

in order to allow for the construction of an effective op- 
timized boson basis (OBB) on each site fc, given by 

l^fe) = E ^s.„. \nk) (nfe = 0, . . . , dopt - 1) . (S8) 

The VMPS ansatz (S6) for the ground state of the Wil- 
son chain is completely analogous to standard finite- 
size DMRG^, and the use of an optimized local basis 
[Eqs. (S7) and (S8)] was pioneered in Ref. 8, finding sub- 
sequent applications in, for example, Refs. 9,10. The 
variational optimization of the resulting MPS with re- 
spect to the Hamiltonian given in the main paper, 
depicted in Fig. 1 there, proceeds by iteratively updat- 
ing the A- and ^-matrices through a series of sweeps 
through the chain. Given the directed structure of the 
Wilson chain from large energy scales (left side of the 
MPS) to small energy scales (right side of the MPS), 
similar to the NRG, variational energy-lowering updates 
are performed only when sweeping from left to right. In 
contrast, during the reverse sweep from right to left the 
physical state (and its energy expectation value) is left 
unchanged. Nevertheless, during the reverse sweep the 
>l-matrices are recast into a right-orthonormalized form, 
to ensure that the right low-energy part of the Wilson 
chain is described in terms of properly orthonormalized 
ejfective basis sets. 



To get started, the MPS is always initialized randomly, 
followed by proper orthonormalization. To update the 
coefficient spaces related to a given site k of the Wilson 
chain, we perform the following steps, depicted schemat- 
ically in Fig. S2 (we suppress the index k below): 

(a) At a given site, the starting point is a well-defined 
local setting in terms of orthonormal basis sets |a), 
1/3), and \fi) of the ^-tensor for an effective left, 
right, and full local state space, respectively, as de- 
picted in Fig. S2(a). The current approximation 
to the overall ground-state wave function therefore 
has its focus on the A-tensor of site n, and, setting 
Aaiin = {A[n\)ai3, Can be written as 

\,p) = Y,A^pnWm\h). (S9) 

(b) Downward-orthonormalization of A-tensor to move 
the focus to the V-matrix: combine the left state 
space a with the right state space j3 into one index, 
and use the effective local state space n as another 
index. Then singular value decomposition (SVD) 
of the resulting matrix, 

Aafjn = = ^(Q^)p'^p-^pn J (SIO) 

P 

generates a new orthonormal tensor A' that de- 
scribes a combined description of the product space 
(a, /?), such that 

(a/3) 



3 



The wave function can now be represented as 



/ ^ ^pBph I. 
ph 



In). 



(S12) 



This description of the rest of the system in terms 
of an orthonormal effective basis \p) is used to rep- 
resent operators that connect the rest of the system 
with the current local state space. Being orthonor- 
mal, it also makes the numerics simple and stable. 

(c) By contracting the remaining box B from the last 
step onto the ^-matrix, this shifts the focus down 
to the y-matrix 



(S13) 



which thus has been altered. Having shifted the 
focus, the ground state is now represented as 



T.^pn\p)\r 



(S14) 



where the F-matrices connect the orthonormal 
state spaces p and with the local boson space n. So 
far, the transformation of 1-0) has been exact with 
Eq. (S14) describing the same state as Eq. (S9). 
Transform the Hamiltonian and operators into the 
basis \p) I n) and solve the eigenvalue problem 



E 

(pn) 



H(p'n'){pn)^{pn) — ^g^(p'n')i 



(S15) 



for the ground state of the system. Thus the result- 
ing ^-matrix represents the updated ground state 
via Eq. (S14). 

(d) Upward-orthonormalization of the V-matrix: the 
singular value decomposition of V, 



(S16) 



provides a new effective description of the local bo- 
son space, such that the transformation matrix V 
from the original boson basis to the OBB is orthog- 
onal: 



E 



V -V*~, 



(S17) 



Note that the singular values Sn indicate the rel- 
ative importance of the optimal boson bases. We 
will come back to this point at the end of this sec- 
tion. At this step the wave function is 



IV') = ^SnCnp\n)\p). 



(S18) 



(e) By contracting the remaining box SnC from the 
previous step onto the A'-matrix, similar with step 
(c), the focus can be shifted to the A-matrix again: 



l^p^n^hp- 



(S19) 



The wave function is thus reexpressed in the same 
form as Eq. (S12). Now transform the local oper- 
ators to the OBB \n) using V , and optimize the 
AajSn matrix in the same way as done using the 
traditional VMPS method. 

(f) Combine the left and local indices and perform a 
singular value decomposition of the Aa/sn matrix: 



AaPh 



Ho 



9/9- 



(S20) 



The resulting tensor Aaqn is left orthogonal: 

^'^AaqiiA'^^if^—Sqq'. (S21) 



np 



Contracting the remaining rqFqp to the A tensor on 
the right side completes the update of the current 
site. 

The OBB method enables us to increase the number dk 
of local states that can be kept in the original basis from 
a few dozen to dk ^ 10'*. (In the next section we shall 
show that by implementing explicit oscillator shifts, the 
effective number of local boson states that are accounted 
for in the unshifted basis can be increased to more than 
IQio.) 

The two adjustable VMPS parameters are the dimen- 
sion D of the VMPS matrices (corresponding to the num- 
ber of DMRG states kept) and the dimension dopt of the 
optimal boson basis. To exemplify the influence of D 
and dopt on physical quantities. Fig. S3 plots the magne- 
tization (ax) for SBMl as a function of a for different D 
and dopt- Clearly, ac is already well converged through- 
out. In practice, we chose D and dopt large enough to 
ensure that all singular values [sn and in Eqs. (S16) 
and (S20)] larger than 10~^ were retained throughout the 
entire Wilson chain, except possibly at its very end. We 
will explain this in more detail in the next section, when 
discussing Fig. S7 below. 

For SBM2 with two bosonic baths, on the other hand, 
the combination of two boson sites into one supersite 
requires numerical parameters such as D and dopt to be 
set to larger values than for SBMl. Nevertheless, we 
find that the number of kept states needed to ensure an 
accuracy comparable to that of SBMl is smaller than 
the 

or f^opt that might have been naively expected 
from the fact that the local state space now has a direct 
product structure. 
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FIG. S3: Convergence check for the VMPS parameters D (a) 
and dopt (b), for SBMl. The panels show {px) as function 
of for different choices of D and dopt, respectively. Since 
critical exponents are obtained from fits to the linear parts of 
such curves, the resulting exponents are evidently not sensi- 
tive to D and dopt. As in Fig. SI, the dashed lines show the 
mean-field power law with /3mf = 1/2. 



III. OBB WITH EXPLICIT SHIFTS 

The OBB scheme described in the preceding section al- 
lows us to easily perform calculations on a desktop com- 
puter using local boson bases of dimension < 10''. 
This can be increased by at least 6 more orders of mag- 
nitude by analytically incorporating explicit shifts during 
the construction of OBB. The idea is to explicitly shift 
the harmonic oscillator coordinates ifc by their equilib- 
rium expectation values (5:^) (similar in spirit to the pro- 
cedure used in Ref. 11), such that the OBB can be used 
to capture the quantum fluctuations of the shifted coor- 
dinate 



{Xk) 



(S22) 



We now describe explicitly how this is done. 

We begin by noting that a shift corresponds to a uni- 



FIG. S4: Convergence check for the VMPS parameters D 
(a) and dopt (b), for SBM2. The panels show the transverse- 
field response, \{o'z}\ vs. in the CR phase, with a robust 
sublinear power law. The deviations at smallest fields arise 
from small symmetry-breaking effects in the numerics which 
cause a linear response akin to an ordered state. 



tary transformation 



UiSk)^e%^^-'^^ 



such that, 



= U\Sk)bkUiSk) = bk + 



Sk^ 
V2 



(S23) 



(S24) 



on bk (similarly for b\). Thus the harmonic oscillator 
displacement Xk = ^(^fc + ^l-) is shifted to x'j, — Xk + Sk, 
and the local boson number operator to 

h'^-b'fb'^ =nk + SkXk + §. (S25) 

The shift can be implemented on the Hamiltonian level, 
by replacing the original Wilson-chain Hamiltonian by 
the shifted Hamiltonian 



y^^'ibM,bi^ = nl^'(b'k^,m 



(S26) 




boson site k 



FIG. S5: Effective boson occupation number {fik) in the 
original, unshifted Wilson-chain boson basis, as function of 
site index k along the Wilson chain, calculated for the lo- 
calized phase of SBMl. Dashed line indicates the relation 
{hk) oc A*'^"^' (see text). The finite size effect causes VMPS 
result deviate from this exponential relation at the end of the 
Wilson chain. 

The local states jnt) in Eq. (S8) now represent Fock 
states of the shifted oscillators. 

To incorporate explicit shifts within the OBB sweep- 
ing strategy described in the last section, we calculate 
the ground state expectation value (xfe) of the current 
site after step (e). We then set 5k = —{ik) in Eqs- (S24) 
and (S25), thereby ensuring that the shifted coordinate 
x'l. corresponds to Eq. (S22). We subsequently move back 
to step (b) and implement the shift in the Hamiltonian. 
(In practice, it is convenient to preserve the form of the 
Hamiltonian itself, and instead to change the matrix rep- 
resentation of the operators bk, 6^ and fik to implement 
the shift of Eq. (24).) Then we repeat the local update 
loop of the sequence (b), (c), (d), (e) (Fig. S2) until the 
shift converges. Finally, we move to step (f ) and the next 
site. 

The SBB method allows us to reach boson shifts (xfe) 
so large that their description within the unshifted boson 
basis would require local dimensions of order ~ 10^*^, 
while nevertheless keeping the actual number of boson 
states in the shifted basis reasonably small, typically 
dk ^ 10^. A typical result for the resulting boson oc- 
cupation numbers (fik) in the original, unshifted Wilson- 
chain boson basis is shown in Fig. S5, calculated in the 
localized phase. Since the boson shifts for the bosons in 
the original definition of SBMl scale as'^ (xk) — 
with uJk oc A"*^, we expect and indeed find that (fik) in- 
creases exponentially with k, as A*^'^""^. 

An accurate representation of this exponential rise, as 
achieved by our VMPS scheme, is essential for obtaining 
correct results for critical exponents. The detrimental 
effects of Hilbert space truncation are illustrated vividly 
in Fig. S6. It shows (a^) vs. {a — ckc) for SBMl, calcu- 
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FIG. S6: The effect of Hilbert space truncation on the calcu- 
lation of the SBMl critical exponents /3 in (a) and 5 in (b). 
The calculation labeled "d^ = 100, shifted" was performed us- 
ing the OBB method with explicit shifts, without setting any 
maximum value for the size of the shifts, i.e. without restrict- 
ing d^c^. Dashed straight lines correspond to power laws with 
exponents /3nrg = 1-18 and 5nrg = 1.85 obtained by NRG^'^. 
<5nrg agrees with the hyperscaling exponent 5 = (l + s)/(l — s) 
that applies for s > 1/2, see Sec. IV C below. 



lated for several upper limits on the size of the local 
boson basis. While the exponents obtained by NRG^'' 
(indicated by dashed lines) correspond to < 12, the 
curves clearly change strongly as is increased. Indeed, 
fully converged results are obtained only when the local 
boson basis can be taken to have "unlimited" size, as is 
the case for the explicit shifting strategy discussed above. 

With the SBB method, the results for chain length 
L = 60 shown in this paper can be obtained within a 
few hours on a desktop computer. Note, though, that in 
the localized phase the total time needed for the calcula- 
tion increases exponentially with L. The reason is that 
the converged value for the effective shift (x^) increases 
exponentially with fc as 

(ifc) ^A'=(i-^)/2, (S27) 



6 



as explained above. However, the "sweeping step size" 
for {xk), i.e. the change in this quantity from one sweep 
to the next, is Umited, in effect, by the dimension dk of 
the shifted local boson basis. Consequently, the num- 
ber of sweeps needed to achieve convergence for (a;^.) in- 
creases exponentially with k (and making an informed 
initial guess for the requisite shift {xk) does not really 
help in speeding up its accurate determination). 

For a given set of convergence criteria, the exponen- 
tial growth in the shifts {xk) is accompanied by a similar 
growth in the absolute errors in their numerical deter- 
mination. The consequences of this for SBMl can be 
seen in Figs. S7(a) and S7(b), which show examples of 
the singular values [Eq. (S20)] and Sh [Eq. (S16)], re- 
spectively, as functions of the Wilson chain k: when k be- 
comes large, the the lower end of the singular value shows 
an increasing amount of scatter. This implies that these 
singular values are not yet properly converged, which 
is directly correlated to the uncertainties in the oscilla- 
tor shifts. Better convergence can be achieved by using 
stricter convergence criteria, as illustrated by Figs. S7(c) 
and S7(d), but only at a considerable increase in com- 
putation time, essentially using up to several hundreds 
of sweeps. We have thus adopted a compromise between 
accuracy and computation time: we chose D, d^pt and 
the shift convergence criteria such that throughout the 
entire Wilson chain, all discarded singular values were 
smaller than 10~^, except possibly for an increase in this 
tolerance at the very end of the chain, of the type seen 
in Figs. S7(a) and S7(b). We have checked explicitly 
that not-optimally-converged shifts and singular values 
towards the end of the Wilson chain do not noticably af- 
fect the resulting physical quantities of interest, i.e. that 
these are already well converged nevertheless. 

Interestingly, we have found that for SBM2 it is easier 
to obtain a well-converged shift in the localized phase 
than for SBMl, because the competition between the two 
chains causes the increase in oscillator shifts near the 
end of the Wilson chain to be much smaller for SBM2 
than SBMl, as illustrated in Fig. S8. Near the chain's 
beginning, in contrast, the singular values for SBM2 were 
found to be bigger than for SBMl, but we nevertheless 
ensured throughout that only singular values below the 
10~^ tolerance were discarded. 



IV. ADDITIONAL RESULTS FOR SBMl 
A. Determination of phase boundary 

We have explored a number of different methods for de- 
termining the phase boundary between the localized and 
delocalized phases, all of which yield essentially equiva- 
lent results. (Fortunately, no oscillator-shift-related con- 
vergence problems occur in this context, since the oscil- 
lator shifts are essentially zero at the phase boundary.) 

1. "Best power law". As mentioned in the main pa- 
per, our "standard method" for determining Uc has been 
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FIG. S7: Panels (a,c) and (b,d) show examples of the sin- 
gular values Vq [Eq. (S20)] and sa [Eq. (S16)], respectively, 
as functions of the Wilson chain k, calculated in the local- 
ized phase of SBMl. Panels (c,d) show the same quantities 
as (a,b), but calculated using stricter convergence criteria re- 
quiring more sweeps, resulting in better-converged singular 
values [see discussion after Eq. (S27)]. The maximum value 
of {hk) in this example is of order 10*, corresponding to a 
maximum shift (Sk) is of order 10''. The diagonal dashed 
lines show the fc-dependence of the shifts 5^, multiplied by 
a constant prefactor that was chosen by hand in such a way 
that the dashed lines he near the onset of noisiness at the 
lower end of the singular value spectra. This reveals the di- 
rect correlation between the exponential increase of 5| , which 
directly enters the Hamiltonian, and the noise in the smallest 
singular values. 
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FIG. S8: Panels (a,b) show examples of the singular values 
rq [Eq. (S20)] and sn [Eq. (S16)], respectively, as functions 
of the Wilson chain k, calculated in the localized regime of 
SBM2. For all the SBM2 parameters studied in this paper, the 
maximum effective boson occupation number (fik) is relatively 
much smaller than those we found in SBMl {{fik} < 10''), 
which results a faster convergence of shift {Sk}- In this typical 
example for SBM2, we get converged singular values already 
with A^swccp — 34. 
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to tune it such that {a^) vs. [a^ — etc) yields the best 
straight hue on a log- log plot. This turned out to be the 
most convenient way of getting accurate critical expo- 
nents. 

2. Energy flow diagrams. In NRG, energy-flow 
diagrams^ can be used to accurately determine the phase 
boundary of a quantum phase transition, since the de- 
generacies of low-lying levels typically differ for the two 
phases to be distinguished. It is possible to gener- 
ate such energy-flow diagrams also within the VMPS 
approach^^, by calculating the eigenstates of the left 
block's Hamiltonian at each site when sweeping from left 
(large energy scales) to right (small energy scales), and 
then appropriately rescaling and shifting the resulting 
eigenenergies^ . Details of this procedure and explicit ex- 
amples of flow diagrams obtained in this manner will be 
presented elsewhere^. 

Energy-flow diagrams, however, are sensitive to rele- 
vant perturbations in terms of small numerical inaccura- 
cies at large energies, i.e. early Wilson shells. This can 
lead to the artificial breaking of symmetries. For exam- 
ple, the Hamiltonian of SBMl at = Q commutes with 
the parity operator 



P = cr^e' 



(S28) 



where N — h\hj. counts the total number of bosons 
on the entire Wilson chain. The corresponding parity 
symmetry guarantees that the ground state is two-fold 
degenerate. However, this degeneracy is typically broken 
by numerical inaccuracies, unless the symmetry is explic- 
itly implemented in the numerical code. We have done 
so, and flow-diagrams resulting from a parity-symmetric 
version of our code yield Uc values in agreement with the 
"standard method" used in the main paper. 

3. Diverging boson number. Another procedure for 
determining ac is to plot the boson occupation number 
for all sites of the Wilson chain: a divergence of (n^) with 
Wilson-chain index fc, as seen in Fig. S5, is a signature 
of the localized phase. 



B. Critical exponents 

Fig. S9 shows some typical data sets used to extract the 
critical exponents /? and 5 for s < 1 shown in Figs. lc,d 
of the main text. This allows to assess the accuracy of 
the VMPS method in the quantum critical regime. 

To properly describe the critical behavior, the Wilson 
chain must be long enough to resolve energy scales down 
to the scale T* oc |a — ad"^, which bounds the quan- 
tum critical regime^; here v is the correlation- length ex- 
ponent. Now, the lowest energy scale accessible for a 
bosonic Wilson chain of length L, i.e. L — 1 boson sites, 
is A~^. Thus, to determine ac with an accuracy of, say, 
10"'*, we need A^^^-^) < cx 10"'^'', implying that the 
requisite chain length scales as 



For SBMl, the correlation-length exponent becomes 
large for small s, see Fig. 5a of Ref. 2, and hence the 
requisite chain length increases with decreasing s, too. 
Together with Eq. (S27), according to which the largest 
shifts at the end of the chain scale as A^'^^"'''/^, this im- 
plies that in the localized phase, the sweep time needed to 
reach convergence increases exponentially as s decreases 
below 1/2. 

For the above reasons, the data for s = 0.2 is clearly 
less accurate than for s > 0.3, as reflected by the er- 
ror bars shown in Fig. 2c, d of the main text. The two 
extremal values that define the indicated error bars cor- 
respond to the two values of the exponent obtained by 
using only the upper or lower half of the full fitting inter- 
val that is indicated by vertical marks; the uncertainties 
from the straight power-law fit are somewhat smaller. 

For s = 1/2, the transition is at its upper critical di- 
mension and logarithmic corrections to the leading power 
laws are expected. As quantifying these corrections from 
the VMPS results is difficult, we have restricted ourselves 
to fits to effective power laws, which then naturally re- 
sult in exponents slightly deviating from the mean-field 
values. 



C. Hyperscaling 

For SBMl, a scaling ansatz for the singular part of the 
free energy can be used^'' to derive a hyperscaling relation 
between the critical exponents 5 and a;, namely (5 = (1 + 
x)/(l — x). Here, x describes the divergence of the static 
susceptibility as T — ?> at criticality, x{T) oc where 
X = d{(Jx) /dh^. Hyperscaling also implies x = y, where 
y characterizes the divergence of the zero-temperature 
dynamic susceptibility, xi'^) ^ uj~'^. Furthermore, y = s 
is an exact result for the critical long-range Ising chain 
at aU s (Refs. 14,15). This finally yields 



1 



1 



(S30) 



In A 



(S29) 



under the condition that hyperscaling is fulfilled, i.e., 
that no dangerously irrelevant operators spoils the naive 
scaling hypothesis. For the Ising chain, this applies be- 
low the upper critical dimension, i.e., for s > 1/2, while 
a dangerously irrelevant operator appears for s < 1/2. 
Consistent with the expectations from QCC, our (5-values 
follow the hyperscaling prediction (S30) for s > 1/2 (see 
Fig. 2d of main text). 



V. ADDITIONAL RESULTS FOR SBM2 

A. Determination of phase boundaries 

Below we describe several approaches that we have 
found useful for determining the various phase bound- 
aries of SBM2. For convencience, these boundaries are 
illustrated schematically in Fig. SIO. 
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FIG. S9: The left and right columns show some of the VMPS 
results used to determine the SBMl exponents /9 and 5 that 
are shown in Fig. Ic) and Id) of the main text, respectively. 
We used = 0.1, A = 2, L = 60, D = 40, dopt = 12. The 
main text (Eq. (S29)) explains why the range of pure power- 
law behavior is smaller for s = 0.2 than for larger s- values. 

1. Order parameter. Measurements of the order pa- 
rameter {(Jx,y) of the localized phase can be used to de- 
termine the critical field /i" that defines the location of 
the L-D boundary at finite h^. This is demonstrated in 
Fig. Sll for s = 0.6 and a = 0.1. The method of "best 
power laws" is then suitable to obtain accurate values for 




FIG. SIO: Schematic sketch of the SBM2 phase diagram as 
function of s, a and h^. The localized phase (L) lies below 
the curved surface representing the critical field hl{s,a), the 
delocalized phase (D) above and to its right. (The curved sur- 
face was not calculated, but is an artist's impression.) The 
cross-hatched area in the hz = plane represents the critical 
phase (CR). Its phase boundary with the localized phase (L) 
is marked by the thick red line in the hz = plane (corre- 
sponding to the red line in main text. Fig. 3a), consisting of a 
straight line between the points (s*,0) and (s*, q:c(s*)), and 
a curved portion representing the line ac(s), for s* < s < 1. 
The thick blue line at s = 1 is the CR-F phase boundary to 
the free phase (F). The thick black line at a = hz = likewise 
represents a free spin (F). 

/i^ as well as corresponding critical exponents - we leave 
this for a future study. 

Notably, the inherent instability of the CR phase, com- 
bined with numerical errors of VMPS, renders difficult 
the accurate determination of the CR-L phase bound- 
ary at /i = 0. In particular, a direct observation of the 
order parameter as function of a or s leads to sizeable 
uncertainties. 

2. Transverse-field response. The zero-field phase 
boundary CR-L is most accurately extracted via the re- 
sponse to a small transverse field, {az){hz). The stable L 
phase responds linearly, while the CR phase responds 
with a non-trivial sublinear power-law, see Eq. (S32) 
below - those can be easily distinguished, as shown in 
Fig. 3b of the main paper. Indeed, by studying the 
transverse-field response for small a and 0.7 < s < 0.8, 
we are able to determine the value of the universal critical 
"dimension" s* to be 0.75 ± 0.01, as shown in Fig. S12. 

The accurate determination of s* requires the numeri- 
cal results to be reliable down to transverse fields as small 
as hz ~ 10"'^'*. Such precision is achievable, in principle, 
by using sufficiently large D and dopt, as follows from 
Fig. S4. However, such a brute force approach is com- 
putationally demanding. A more efficient strategy is to 
exploit parity symmetry: For h^ ~ hy — the SBM2 
Hamiltonian commutes with the parity operator 

Pz = aze'^^, (S31) 

where N = Y.,^ blJ)^^ + Y.k bLky counts the total num- 
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FIG. Sll: (a) Order parameter |(cr2:,j,)| for the L-D transi- 
tion, driven to zero by increasing the transverse field past 
the critical value h% (indicated by black dashed line), (b) 
Correspondingly, the transverse-field response of \ {(Jz)\ shows 
a slight kink, indicating the higher-order singularity expected 
at the L-D transition. As in Fig. 3 of the main paper, the 
linear response of the L phase is clearly visible. We observe a 
small breaking of the XY symmetry of SBM2 in the L phase 
due to numerical errors, in that the order parameter in panel 
(a) prefers configurations with \{(Jx)\ = {{o'y)], instead of ex- 
ploring the full rotational symmetry. 



ber of bosons on the entire Wilson chain. Thus, the 
ground state has a parity degeneracy. In fact, this de- 
generacy is the main source of numerical error for small 
hz when using a VMPS code that does not account for 
parity symmetry, since then the degeneracy is lifted by 
numerical noise (similarly to the SBMl case). By imple- 
menting this parity symmetry explicitly in the code, we 
were able to achieve the required precision (much better 
than in Fig. S4), while using choices for D and dopt that 
were not unreasonably large, as shown in Fig. S12. 

3. Diverging boson number and flow diagram. We have 
also explored alternative ways for determening the L-CR 
phase boundary, based on monitoring the the divergence 
of the boson number per site at the end of the chain. 




FIG. S12: Transverse-field response for SBM2, 1((t^)| vs. h^, 
used to determine the L-CR phase boundary. The fitting ex- 
ponent 5' is defined in Eq. (S33), and the fitting range is 
G [10~^^, 10~*]. From the transition between a linear re- 
sponse for s < s* and a sublinear power law for s > s*, we 
estimate s* to lie between 0.74 and 0.76. To reduce numerical 
noise, the curves were calculated using a VMPS ground state 
that also was an odd eigenstate of the parity operator Pz of 
Eq. (S31). 



or analyzing NRG-like flow diagrams (as discussed for 
SBMl). We have found these methods to be computa- 
tionally much cheaper than studying the transverse-field 
response, while yielding results of comparable accuracy 
for the L-CR phase boundary. 



B. Properties of the critical phase 

The CR phase corresponds to a partially screened (or 
fractional) spin, with non-trivial power-law autocorrela- 
tions of the components of a. In ground-state calcula- 
tions, those can be probed by measuring the response 
to an applied field: The linear-response susceptibility at 
T = is infinite, and the non-linear response is of power- 
law character, 



(S32) 



with 6,6' > 1. A standard renormalized perturbation 
expansion around the free-spin fixed point results in^' 



1/6 = 



1 - s 



-0{[l 



1/6' - l-.s + 0([l-s]2). 



(S33) 



In Figs. S13 we show numerical VMPS data for the 
non-linear response for several sets of parameters inside 
the CR phase. We indeed find the expected power laws 
over a sizeable range of fields; the power laws are cut off 
at small fields by the influence of numerical errors, like 
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FIG. S13: Non-linear response inside the CR phase. We find power law behavior as expected from Eq. (S32). For each s, we 
chose an a- value close to the a* fixed-point value, to ensure that the asymptotic power law is reached quickly (i.e. without the 
need for to become extremely small). Using the fitting range between 10~® and 10"^, we find values for the exponents (5flt 
and (5flt in excellent agreement with the values 5pcrt and S'^crt expected from Eq. (S33), as indicated in the legends. 



the breaking of the XY symmetry of the model or the 
generation of a small transverse field. From the data we 
can extract exponents 5 and 5' as function of s. Compar- 
ing the resulting exponents to the perturbative prediction 
(S33), we find that the agreement is excellent, considering 
that (i) for s very close to unity the asymptotic regime is 
difficult to reach, due to the slow flow of a, and (ii) for 
larger (1 — s) the second-order terms in Eq. (S33) (i.e. 
two-loop corrections) will become important. 



C. RG flow near s* 

Fig. 4 of the main text shows the RG flow of SBM2 
in the two cases (a) s* < s < 1 and (b) < s < s*. In 
case (a) the CR phase is stable for small a and ft, = 0, 
whereas this phase is absent in case (b) . In the following 
we briefly discuss the evolution of the RG flow from case 
(a) to case (b), as this is related to a discontinuous jump 
of the phase boundary upon varying s. 

As described in the main text, the fixed points CR 
and QCl of case (a) approach each other upon lowering 
s, such that they meet at s = s* and both disappear for 
s < s* . This follows from the absence of the CR phase in 
case (b) and the fact that both are intermediate-coupling 



fixed points for s < 1 (i.e. they have to meet at finite a). 
The merger of the two fixed points also implies that the 
flow lines to the left and right of CR in Fig. 4a merge. As 
a result, the RG flow line which lead from F via CR to 
D in case (a) now becomes a flow line from F to QC2 in 
case (b). For the phase boundary of the L phase (thick 
line in Fig. 4) this means that its starting point on the 
hz = axis jumps from a finite value (corresponding to 
QCl) in (a) to zero in (b) once s is lowered past s* . 

We note that the merger and disappearance of two 
fixed points upon variation of a "dimension" d is not un- 
usual. For instance, upon approaching the lower critical 
dimension d~ of a magnet from above, the critical fixed 
point typically approaches the trivial fixed point describ- 
ing the ordered phase, such that both merge at d = d~ 
and disappear for d < d~ . Here the evolution of the phase 
boundary is continuous. What is unusual about SBM2 
is that two intermediate- coupling fixed points merge and 
disappear, causing the discontinuous behavior. The only 
other example with similar physics we are aware of is in 
the two-channel Anderson/Kondo impurity model with 
power-law density of states cx where a critical fixed 
point merges with a stable non-Fermi liquid fixed point 
at some critical dimension rmax, with a consequent jump 
in the phase diagram. ^^'^^ 
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